Coherent Backscattering of Ultracold Matter Waves: Momentum Space Signatures 
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Using analytical and numerical methods, it is shown that the momentum distribution of a matter 
wave packet launched in a random potential exhibits a pronounced coherent backscattering (CBS) 
peak. By analyzing the momentum distribution, key transport times can be directly measured. The 
CBS peak can be used to prove that transport occurs in the phase- coherent regime, and measuring its 
time dependence permits monitoring the transition from classical diffusion to Anderson localization. 

PACS numbers: 03.75.-b, 05.60.Gg, 42.25.Dd, 72.15.Rn 



Disorder has dramatic effects on the quantum trans- 
port of matter. Spatial randomness and phase coherence 
together can completely suppress diffusion, as demon- 
strated by the paradigmatic phenomenon of Anderson 
localization During the past decade, there has been 
a growing body of evidence for three-dimensional (3D) 
localization in random media with different types of non- 
interacting waves: light microwaves [3] and ultra- 
sound 0. This ubiquitous and yet elusive phenomenon 
has sparked considerable interest in the field of ultracold 
atoms [il-Q- Key experimental achievements include ID 
Anderson localization in speckle [8i] and quasi-periodic 
[il potentials, as well as 3D localization in momentum 
space with the kicked rotor [lo| . Recently, 3D Anderson 
localization of noninteracting ultracold fermions and 
bosons [l2| in a laser speckle field was reported. 

To claim Anderson localization, one needs to discrimi- 
nate interference-induced absence of diffusion from classi- 
cal trapping or slow diffusion. This requires evidence for 
phase-coherent transport. Here, the coherent backscat- 
tering (CBS) phenomenon is of key importance because 
it arises by interference of waves in random media and 
measures mesoscopic phase coherence [13]. With classi- 
cal waves, CBS appears as an enhancement of the dif- 
fuse intensity reflected off a disordered medium around 
the backscattering direction, and has been observed in 
numerous experiments involving light [3, [HI , but also 
acoustic [l6| and seismic waves [13 • The interference 
causing CBS is also responsible for weak localization, by 
reducing the diffusion coefficent compared to its phase- 
incoherent, or classical, value [l, IHI. In electronic sys- 
tems, weak localization is invaluable for a careful char- 
acterization of phase coherence [l8|. With cold-atomic 
clouds expanding in random potentials, however, the dif- 
fusion constant extracted from real-space data hardly 
shows clear evidence of localization corrections, because 
the cloud contains many different momenta that combine 
to a rather involved spatial profile [19, 20]. 

With this Rapid Communication, we propose to study 
the dynamics of a matter wave that is launched with 



an initial momentum larger than its momentum spread 
in the bulk of a 2D or 3D random potential. Combin- 
ing a numerical and theoretical analysis, we show that 
the CBS signal can be directly observed in the momen- 
tum distribution and studied as a function of time. Ul- 
tracold atoms are an invaluable asset as they offer the 
unique opportunity to visualize the CBS effect on the mo- 
mentum distribution directly, measured inside the disor- 
dered medium. With this setup, one avoids the boundary 
conditions that severely complicate both theory and ex- 
periments of wave scattering by random media [l3|, [2l| . 
Moreover, the momentum-space analysis gives immediate 
access to key mesoscopic parameters such as scattering 
and transport times. Lastly, we demonstrate that the 
CBS measurement provides precious information on the 
phase coherence of the matter wave, and finally permits 
to monitor the transition from diffusion to localization. 

Let a cloud of atoms with mass m be prepared at 
time t = in the state ^(/c, 0) describing a wave packet 
with mean momentum hko = mvo and small spread 
Ak <C \ko\. This can be achieved by releasing the 
atoms from a shallow trap, and either launching them 
with mean velocity Vq or moving the disorder potential 
with —vo relative to the cloud. We assume negligible 
interaction effects. This is the case in practice for a 
very dilute cloud Q or spin-polarized fermions [ll|. For 
concreteness, as realized by harmonic trapping of non- 
interacting particles, we then take the initial distribution 
Po{k) = |^(/c,0)p to be an isotropic Gaussian, 

po{k) = (2^Ai^-2)^/2 exp [-{k - /Co) V2Afe'] , (1) 

normalized to Trpo = / po(^)d/c/(27r)^ = 1. From time 
t = onwards, the matter wave then evolves according 
to the Schrodinger equation with single-particle Hamil- 
tonian H = /2m + V{r). A well-controlled random 
potential is provided by laser speckle [22|]. Without 
loss of generality, V{r) = 0, where the overbar denotes 
the ensemble average over disorder realizations. The 
random potential is then characterized by its variance 
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FIG. 1. (Color online) Momentum distribution 'p{kx,ky,t) of 
a matter wave packet launched with initial momentum ko = 
(A;o,0) in a 2D random potential with correlation length 
averaged over 960 disorder realizations. The time unit is = 
m(^^/h. (a) t = 10t(^: elastic scattering depletes the initial 
wave packet, centered at ko (peak values not shown), and 
populates the disorder-broadened energy shell along the circle 
\k\ = \ko\, while the CBS peak emerges at —ko. (b) t = 18r^: 
The CBS peak is now the dominant feature, proving phase- 
coherent multiple scattering. 



V(r)^ = and spatial correlation length (. This length 
defines a correlation time = m('^/h and a correlation 
energy E^ = h'^/{mC^) ||. 

First, we study the dynamics of the matter wave ex- 
panding in a repulsive 2D speckle potential by solv- 
ing numerically the Schrodinger equation for a potential 
strength V = E(^^ and initial condition /cqC = 2, Ak = 
0.01A:o/V2. For ^'^Rb atoms in a correlated potential with 
C = 1 /im, the absolute time scale is = 1.365 ms. Fig.[T] 
shows the numerically computed, ensemble-averaged mo- 
mentum distribution 'p{kx,ky,t) at two different times. 

At short times [Fig.IlJa), t = 10r<^], one sees a very nar- 
row peak at /co, a broad, ring-shaped anisotropic back- 
ground and a rather smooth peak at —ko. The forward 
peak is the remainder of the initial momentum distribu- 
tion, which is depleted because atoms are scattered out 
of the initial mode at a rate given by the elastic scat- 
tering time Ts The latter can be extracted from 

the early-time decay p(/co, ^) ^ |^(^o,^)| exp(— t/rg). 
We find = 1.43r^, with corresponding mean free path 
is = vo'Ts such that k^is = 5.72 for the parameters used. 
Weak-disorder perturbation theory [6] predicts too low a 
value {kois = 2.32), as known for rather strong, spatially 
correlated disorder This shows how the early-time 
momentum distribution can be used to measure the key 
parameter Tg, even in the strong disorder regime where 
precise analytical predictions are not available [23]. 
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FIG. 2. (Color online) CBS peak contrast C (blue circles) and 
angular width (red squares) as obtained from the numer- 
ics after averaging over 960 disorder configurations and over 
a time window of lOr^. Solid curves: theoretical predictions 
(|1QP and respectively. Dashed horizontal line: angu- 

lar width A^o = 0. 01/^2 of the initial momentum distribu- 
tion, asymptotically reached by the CBS width for times much 
larger than the coherence time ta = {2DAk^)~^ = 216t^. 

Atoms scattered out of the initial mode populate all 
other accessible /c-space modes on the energy shell and 
thus appear along the circle \k\ = \ko\ in Fig. [H Due to 
disorder broadening, the energy shell has a finite width, 
of order ij^^ which is larger than the initial width Ak for 
the chosen parameters. After a time of the order of the 
Boltzmann time tb = 8-5r^ (A^o^B = 34) [6], the dynam- 
ics turns from ballistic to diffusive. As the memory about 
the initial direction of propagation gets erased, the dif- 
fusive momentum distribution then becomes isotropic on 
average. More precisely, for t > r^, i.e. when diffusion is 
fully established, we find that the decay of the anisotropic 
Fourier components of the background is well fitted by 
exp(— t/r), where r is the transport time that governs 
the diffusive dynamics. For the present parameters, one 
finds r = 5.8r(^ (corresponding to a transport mean free 
path koi = 23.2). As expected, r is smaller than tb, 
due to weak localization (WL) corrections arising at 
early times and caused by very short CBS loops. In Fig. 
[U^a), the peak at —ko is the incipient CBS signal. 

At longer times [Fig. [IJb), t = 18r^], the initial 
state is totally depleted, the diffusive background is fully 
isotropic, and the CBS peak is the dominant feature. 
Both its contrast C, defined as the height above the dif- 
fusive background, and its angular width A6 slowly de- 
crease with time, as shown in Fig. [2j 

For a quantitative understanding of these observations, 
we now turn to the analytical description of matter wave 
dynamics in dimension d = 2,3. The ensemble- averaged 
momentum distribution p{k\t) at time t is given by [6[ 

p{k\t) = I I ^$^^,^(0,t)po(/c), (2) 

where the intensity propagation kernel ^hk'E projects 
the initial momentum k on the energy shell describes 
the ensuing unitary dynamics generated by the Hamil- 
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tonian iJ, and projects back onto the final momentum 
k' . For long enough times t ^ tb, but well before 
the onset of Anderson localization, the atomic dynamics 
is diffusive with an energy-dependent diffusion constant 
D{E) = 2Et / imd) that incorporates the short-range WL 
corrections. The intensity propagation kernel then takes 
the form 



A{k,E)A{k',E) 
2tiu{E) 



eM-D{E)qH]. (3) 



The spectral function A{k,E) = 27t {k\ 6{E - H) \k) is 
the average probability density that a plane-wave state 
\k) has energy E. It also determines the average den- 
sity of states u{E) = J A{k, E)dk/ {27rY^^ . Using (O at 
momentum transfer gr = in Eq. (j2j) results in a time- 
independent isotropic diffusive background: 



dk f dE A{k,E)A{k\E) 



(27r)^ J 27r 27ru{E) 



Po{k). (4) 



This contribution to the momentum distribution de- 
scribes scattering processes that do not rely on long-range 
phase coherence; the label "L" refers to the ladder topol- 
ogy of the underlying Feynman diagrams, Fig. Ha) [131. 
Phase-coherent multiple scattering leads to an additional 
contribution, given by the maximally crossed diagrams 
shown in Fig. ^h). These diagrams describe the inter- 
ference of amplitudes that propagate along the same set 
of scatterers, but in opposite directions, and give rise to 
the CBS effect. Its contribution 'p^^\k\t) is obtained 
from Eq. (j2j) after the substitution q = k^k' in Eq. (j3j). 



(C) 

kk'E\ 



(0,i) = *Lfc^'£;(fe + fc',t). 



(5) 



Because this coherent kernel peaks at k' = — /c, the CBS 
signal will appear at —k^ if the initial distribution is cen- 
tered around /cq, as borne out beautifully by Fig. [TJ For 
an initial plane- wave state (A/c = 0), the CBS contribu- 
tion at exact backscattering is at all times exactly equal 
to the background level, and thus enhances the diffusive 
density at that point by a factor of two. This has been 
observed in optical and acoustical experiments [3, IHI . 

If, however, the initial state has a finite momentum 
spread, the CBS signal is reduced because Eq. ([2]) con- 
volves the plane- wave kernel (|5]) with the initial momen- 
tum distribution. In order to be able to detect other 
dephasing processes, we therefore need to analyze in de- 
tail the CBS signal for matter wave packets with finite 
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FIG. 3. Ladder (a) and maximally crossed (b) multiple 
scattering Feynman diagrams, giving rise to the diffusive and 
CBS contributions, Eq. Q and Eq. (O, respectively. 



momentum spread. From here on, we assume 
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The first inequality implies the weak-disorder condition 
/co^s ^ 1 [13], whereas the last one requires the atom 
coherence length A/c~^ to span many mean free paths, a 
necessary condition to observe the CBS peak (see below). 

Let us first discuss the diffusive contribution (j4]). The 
regime of interest (|6]) implies A/c <C ^7^, i.e. an initial 
distribution p^ik) much narrower than the spectral den- 
sity. This is apparent in Fig.llja) where the initial peak is 
much narrower than the diffusive background. Therefore, 
the result of their convolution in (|4]) can be approximated 
by A{ko^E) « A^{Eo - E) = hj/[{Eo - Ef^n^j^A], 
where 7 = r'^ and Eq = Eu^ = h^kl/{2m) The 
remaining integration over E gives 



,(L) 



A2^{Ek' — Eq) 



(7) 



The background reaches its maximum p^^^(/co) = 
Ts/[Khv{EQ)\ on shell. Spectral broadening results in a 
relative half-maximum width of Hj/Eq = 2/(/co^s). 

Consider now the CBS contribution, Eq. ([2]) with ([5]). 
Since D{E) is a smooth function of E on the scale of ^17, 
one can approximate 



— (t) 



2nv{Ek) 



ft 



(8) 



The exponent on the right-hand side introduces the 
CBS /c-space width A/cc = l/y/2Dt^ which is the in- 
verse of the diff usive sp read in real space. Written 
as A/cc = ^rd/ (2t), this width is seen to be al- 
ways smaller than the spectral width in the diffu- 
sive regime t ^ tb- It is therefore never permissible to 
approximate the spectral function {Ek — Ek' ) by its 
clean limit 2iT5{Ek — E^') when calculating the CBS peak 
resulting from a wave packet of finite width. 

Within the regime delimited by inequalities (|6]), the 
CBS contribution at k' = —ko-\-q around backscattering 
then follows directly by Gaussian integration. 



p(C)(q_fco,i)=^li^)(feo)x 



exp[-gV2Ag2(t)] 



(9) 



Here, ta = {2DAk'^) ^ is the time it takes the matter 
wave to spread diffusively over the coherence length. The 
CBS contrast C{t) = p^^\-kQ,t)/p^^\ko) is found to 
decrease like 



C(t) = (l + t/rA) 



-d/2 



(10) 



The CBS signal can only be observed if this contrast is 
not too small at time r. The corresponding condition 
r A ^ T requires the initial wave packet to be coherent 
over a distance Ak~^ ^, as stated in (|6]). As function 
of gr, the diffusive CBS signal (|9]) is an isotropic Gaussian. 
Its variance Aq'^{t) = + {2Dt)-^ = Ak'^{l + r/^/t) 



4 



is the sum of the variances of the initial momentum dis- 
tribution ([1]) and the interference kernel ([8]). Expressed 
as a function of the angle 6 away from backscattering at 
fixed radius k' = ko^ the CBS peak width decreases as 

AO{t) = AOo^/lTV^t, (11) 

where A^o = Ak/ko is the angular width of the initial 
state in momentum space. In Fig. [21 the predictions (p!Q|) 
and (pTj) for d = 2 are plotted together with the numerical 
data (ta = 216r(^). The agreement is excellent in the 
diffusive regime, validating our analytical description of 
matter wave CBS in bulk disorder. 

We now discuss the behavior of the CBS peak at the 
Anderson localizatio n ons et. In 2D, this occurs when 
the diffusive spread y2Dt reaches the localization length 
^ = exp(7r/co^B/2), i.e. around the time rioc = i^/D. 
This regime is not visible in Fig. [2] since rioc = lO^^r,^ 
for the parameters used. Applying the self-consistent 
theory of localization [^J], we predict two possible sce- 
narios. If the initial coherence length is much smaller 
than the localization length, (A/c)~^ <C ^, the CBS 
peak width saturates at A6>o already around the time 
Ta <C Tioc- The contrast continues to decrease like 
{TA/ty^^ until it reaches Cioc = {ta/tiocY^'^ <C 1. If, 
on the contrary, the coherence length exceeds the local- 
ization length, (A/c)~^ ^ ^, then the CBS interference 
becomes sensitive to the localization of the wave ampli- 
tudes around the time rioc ^ •, and the width freezes at 
^^loc = (^oO~"^5 contrast, Eq. (p!Q|) . of order unity. 
Observing this effect in an actual experiment would pro- 
vide compelling evidence for Anderson localization, and 
the CBS peak could serve as an independent measure for 
the localization length. Supplementary evidence would 



be gained by perturbing the phase coherence deliberately, 
for example by varying the speckle potential over the 
course of time, in order to dephase the CBS signal and 
suppress localization in a controlled manner. 

Let us conclude by showing that the proposed exper- 
iment is immediately feasable with current techniques. 
Ref.[l2| uses condensed bosons at finite temperature, 
with a momentum spread (A/c)~^ of thermal and con- 
densed components that we estimate at 2.5 /im and 
1.6 /im, respectively. To ensure a good CBS contrast, 
these values should not be smaller than ^, according to 
Eq. dnj. For our choice of /cq = 2/(", we find £ 2.3 /im 
and /cq"^ 0.1 /im for ( ^ 0.2 /im of Ref. [l2|, such that 
k-'^ < £ ^ (Ak)-^. We thus predict that for a slightly 
more coherent cloud or a slightly smaller ko under the 
same conditions, CBS should be clearly observable. Fi- 
nally, we briefly comment on atomic interactions, leaving 
their detailed study for future investigation. Their pri- 
mary effect is to decrease the CBS contrast [2l|. For a 
ratio of order 10 between disorder strength and atomic 
interaction energy, as in Ref. [l2| , preliminary numerical 
calculations show a reduction of C by only 10% at time 
t = 18r^. The CBS effect of matter waves is therefore 
a robust phenomenon that offers a unique opportunity 
to characterize key transport parameters and to monitor 
the transition from diffusion to Anderson localization. 
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